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Abstract 

A useful finite-dimensional matrix representation of the derivative of periodic functions is obtained 
by using some elementary facts of trigonometric interpolation. This N x N matrix becomes a projection of 
the angular derivative into polynomial subspaces of finite dimension and it can be interpreted as a generator 
of discrete rotations associated to the z-component of the projection of the angular momentum operator in 
such subspaces, inheriting thus some properties of the continuum operator. The group associated to these 
discrete rotations is the cyclic group of order N. 

Since the square of the quantum angular momentum L 2 is associated to a partial differential boundary 
value problem in the angular variables 8 and ip whose solution is given in terms of the spherical harmonics, 
we can project such a differential equation to obtain an eigenvalue matrix problem of finite dimension 
by extending to several variables a projection technique for solving numerically two point boundary value 
problems and using the matrix representation of the angular derivative found before. The eigenvalues of 
the matrix representing L 2 are found to have the exact form n(n + 1), counting the degeneracy, and the 
eigenvectors are found to coincide exactly with the corresponding spherical harmonics evaluated at a certain 
set of points. 



PACS: 03.65. Ca, 03.65.Lx, 02.60.Jh, 02.60.Lj 



1. Introduction 

A Galerkin-collocation-type method based in a iV-dimensional matrix representation of ^ obtained 
through Lagrange's interpolation has been used to solve one-dimensional boundary value problems (see Refs. 
[l]-[4]). This technique consists basically in the substitution of the variable x and the derivative 4- by N x N 
matrices X and D in a certain form of the given differential equation. The diagonal matrix X has the N 
different values Xj, j = 1, 2, • • • , JV along the main diagonal, and the matrix D is given by 



(1) D = PDP~ 1 , 

1 



where 



{ N 



(2) Dij 



\- 1 . . 

Pij = p' '(x^Sij , i, j = 1,2,- • -, N, 



. (xi — Xj) 



The symbol J2' appearing in (2) means the sum over I ^ i and the prime on p means differentiation of the 
polynomial 

N 

p( x ) = Y[(x-x k ). 
fc=i 

The N nodes Xj can be chosen by imposing a condition where the coefficients of the differential equation 
and the boundary conditions play the main part. This condition is 

N 

(3) E'r^^-^rT' ; = m,-,at, 



where ■j(x) is a function defined by the boundary conditions and the differential equation (see Refs. [4]-[5]). 

More precisely, the projection scheme for the fc-th derivative of a real function / evaluated at different 
(but otherwise arbitrary) points Xi, i = 1, 2, . . . , N is given by 



N 



(4) f k \x j ) = Y,D k j if{x l )+E h j = l,2,-..,N, 

i=i 

where Ej is the j—th component of the residual vector depending on /, N and k. Let us denote by ttn-i the 
space of polynomials of degree at most N — 1. Thus, if / e 7Tjv-i, it is found that (4) is exact at the nodes, 

1. e., Ej = 0, and therefore, any differential problem closed in 7Ttv-i can be solved also in R N as a matrix 
problem yielding the same solution. An estimate of the error Ej for other kind of functions is given in [4] . 

The formal application of this elementary technique to some multivariate problems is straightforward, 
as shown in the following section (see also [6]). 

2. Multivariate case 

To illustrate how this method should be applied to some boundary value problems in several variables 
we will take first two real variables, x and y. 

Let II = 7TM-1 ® ttjv-1 be the space of bivariate tensor-product polynomials of degree N — 1 in the 
variable x and M — 1 in y. Thus, if / e II, f(x, y) can be written as 

JV-l M-l 

(5) f(x,y) = a JkX j y k - 

j=0 k=0 

Now let us take two sets of different but otherwise arbitrary points {xi,X2, ••• , £jvi } and {yi, t/2, ■ ■ ■ , yjv 2 }, 
on the x and y axes respectively, and let D x and D x be the N x N and M x M matrix representations of 
and constructed according to (1) with such sets. 

By deriving partially (5) n times with respect to x and m times with respect to y, and evaluating the 
result at the cartesian nodes (xj,yk) we obtain the vector f(™<" 1 ) of dimension N = NM whose entries are 
given by 

j(n,m) _ 9 n+m f{x,y) 



dx n dy Ti 
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where the indexes r, j and k are related through 

(6) r = j + (k - 1)N, j = l,2,---,N, fc = l,2,..-,M, 

in such a form that r = 1, 2, • • • , N. It is known that a bivariate interpolation on the grid (xj,yk) is uniquely 
possible in ft (see [7] ) so that / can be taken as a bivariate and sufficiently differentiable real function other 
than a polynomial, and (5) as its corresponding (tensor) Taylor polynomial. By choosing j (the ai-index) to 

run faster than k and using (4) according to the case, we can write down a matrix formula in R N for f ( n ' m ) 
in terms of f = f^ 0,0 ', generalizing (4): 

(7) f (n ' m) = {D y m ® D x n )f+E. 

Here, the Kronecker product of matrices A — (a jfe ) and B — (bjk) of sizes N x N and M x M, respectively, 
is defined by 

/ a\\B ai 2 B ■ ■ ■ cl^nB \ 
a-i\B CL22B ■ ■ ■ a,2, nB 



A®B 



\cijvi-B apf 2 B ■■■ a 7v nB ) 



The components of the residual vector E arc zero if / e ft and an expression for them is given in [8] for 
sufficiently differentiable functions. 

Let l x and l y denote the identity matrices of dimension N and M respectively, and define the N x N 
matrices D^, and T> y by 

(8) T) x =ly®D x , n y = D y ®i x . 
Due to the properties of the Kronecker product, these matrices commute: 

T) x D y = T> y D x = D y ®D x . 

More generally, 

T> x n T> y m = {l y ®D n x ){D™®l x ) = D™®D n x , 
and therefore, (7) takes the form 

(9) f ("»"») = D x "Dj / m f + E, 

indicating that the partial derivatives d/dx and d/dy take in II the tensor-product forms given in (8). 
By using the properties of the tensor product and defining P = P y <x> P x and 

D x = 1„ ® D x , l5 y = D y ® l x , 

where P x , P y , D x and D y have the structure given in (2), it is possible to give the following alternate form 
of (9): 

f(",m) = pD x ™D y m p-lf +E . 

On the other hand, the projection of the coefficient functions of the differential operator can be written 
as diagonal matrices since the the product of a function a(x, y) by the partial derivatives of the unknown 
function evaluated at the nodes (xj,yk) is 

a r f^ m) =a(x j ,y k )f^ m \x j ,y k ) 

(a(xj,yk) should be well defined) and the indexes can be ordered according to (6), producing that in this 
scheme, the coefficient functions can be represented by N x N diagonal matrices whose non-zero elements are 
given by a r — a(xj,yh), where r, j and k are related by (6). Let us denote this generic coefficient matrix by 



A. Thus, the part of the differential operator consisting in the product a(x,y)f^ m ' n \x,y) takes the matrix 
form 

AD/D/f. 

If a(x, y) accept a Taylor expansion (this condition is too restrictive and it can be relaxed, but it is adequate 
for our illustrative purposes), we have that A can be defined by the same function a(x,y) through 

A = a(X,Y), 

where X and Y are the matrices given by 

where the N x N diagonal matrix X has the set of points Xj, j = 1, 2 • • • , N along the main diagonal whereas 
the M points yj lie along the main diagonal of the M x M matrix Y . X and Y represent the variables x 
and y, respectively. 

The generalization to the case of q variables x\, X2, • • • , x q is straightforward. We will choose Nj points 
on the a^-axis, and the projection space as the tensor product of the subspaces of polynomials of degree at 
most Nj — 1 in Xj, i.e., 

q 

II = <g> it Nj -i- 

The nodes will be ordered in such a way that a function f{x\, ■ ■ ■ ,x q ) evaluated at the nodes yields the 
vector f whose entries are 

fr = f{x 1 j 1 ,x 2 j 2 ,---,x q jq ), j k = l,2,---,N k , k = 1,2, • • • ,q, 

and the indexes r and j k are related through 

r = 3i + (.72 - l)JVi + (.73 - l)JViJV 2 ■ • • + (j q - 1)N/N g , j k = 1, 2, • • • , N k , k = 1, 2, • • • , q, 

and TV = rife=i ^fe- The index than runs faster is j\, then j 2 and so on, yielding that r = 1, 2, • • • , N. 
The N x N matrix representation of d/dx k , is now 

D fe = 1, <g> • • • <g> l fe+ i (g) D k (g) l fe _i • • • (g) li, 

where lj is the Nj x Nj identity matrix and D k is a matrix of dimension N k having the structure given by 
(1). Similarly, the representation of the variable x k is 

X fe = l q <g> • • • ® l fe+ i ® X k <g> l fe _i • • • ® li. 

In the next section we adapt this technique to other kind of subspaces and it will be applied to an 
important physical problem in Section 4. 

3. Discrete rotations 

In this section we show that a matrix representation of the angular derivative (the derivative of 
univariate periodic functions), can be related to a generator of discrete rotations (associated to a subgroup 
of the rotation group) and becomes a finite-dimensional matrix representation of the ^-component of the 
angular momentum. 

The notation and language used in this section is that of Quantum Physics, what can be seen as a 
digression in this numerical look-like paper, but we think that the reasons are obvious. 

Let us begin by considering a complete set of quantum states in a space of finite dimension, \<fj), to be 
determined later. Here, ipj indicates the j-th eigenvalue of the operator associated to the spatial observable 
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ip (an angular variable). Now, we ask for the operator that produces a cyclic permutation, within a phase 
shift factor, of the complete set of states • • • , \<Pn), that is, the operator A that yields 



A|^) = 



171 bi), J=N. 

Therefore, the representation of A in the <^-basis is the matrix of elements ((pj\A\(pk) = Ajk given by 



(A jfe ) 



/ 





••• 





e^ 1 \ 


g»72 





••• 













••• 














••• 








V o 





••• 




o / 




A N 


= e< r l, 







(10) 

yielding that 
(11) 

where T — X^jLi 7j> an( l 1 i s tne identity matrix of dimension N. Since the determinant of A is 

detA = (-l) JV+1 e ir , 

A represents a proper rotation and generates a finite subgroup of the rotation group if 

2Ztt, N odd, 

(21 + 1)tt, N even, 
with / integer. In such a case, (11) becomes 

(12) A N 



showing that the matrix set of powers of A is a representation of the cyclic group of order N or 2N for 
N odd or even respectively (being a two-valued representation of our finite group of rotations in the latter 
case) . Therefore, every eigenvalue of any power of A is a rational root of unity. Let us make clear this point 
because it is related with the main result of this section. By fixing the phase shifts jj, we can write (10) 
in more tractable forms, but due to (12) it is necessary to consider the cases of N odd or even separately. 
Thus, we can take the N x N basic circulant permutation matrix 




A,- fc 



for (10) in the odd case, whereas, for N even, 



S jN , for j = l, 

Sj,k+i, for j > 2, 

-dj N , for j = 1, 

8j,k+i, forj>2. 
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It is not difficult to see that the characteristic polynomial of A in the variable A is (— \) N + 1 = for any 
TV, and that the TV-tuple of components e ~ lkaj /\/N, k = 1,2, • • • , TV, is the j-th normalized eigenvector of 
A yielding the eigenvalue Xj = e laj , where 



(13) 
and 



27TJ . T 

a j = — > 3 e In, 



0,±l,±2,---,±n, 



TV = 2n + 1, 



IN 



(14) 



±l/2,±3/2,---,±(2n-l)/2, TV = 2n. 
Being A unitary, it defines an hermitian matrix A through 

A = e iA 



whose eigenvalues are given by (13) (note that the trace of A vanishes) and its eigenvectors, denoted by \aj), 
are those of A. Therefore, the unitary matrix diagonalizing simultaneously A and A is the one with entries 



(15) 



(Vk\aj) 



1 

Vn 6 



-ika i 



fc,j = l,2,---,TV. 



In the usual quantum case A becomes proportional to the angle of rotation times L z . To find a discrete 
analogue relation we need to calculate the elements of A in the y>-basis. This can be made through 



N 



A jk = (<pj\A\<p k ) = ^2(<Pj\ai)ai(ai\ip k ) = ^ ]T le^~ k ^ a '. 

i=i iei N 

Note that the diagonal elements vanish all of them. For N = 2n + 1 we have 

4tt . 2-KUj-k) 



N 



whereas for N = 2n, 



M = -^('-W^'- 1/2)0 -*» 



N 



These sums can be calculated easily [9] and, for iV odd or even, they become 

0, j = k, 

(16) A jk = (<pj\A\<p k ) = { 



i(-iy +k ir/N 



sin 



N 



J ^ k. 



We proceed now to establish a relation between A and a TV-dimensional matrix representation of the 
derivative of trigonometric polynomials. We choose the odd case TV = 2n + 1 first. 



3.1 Odd case 

Let t„ be the space of trigonometric polynomials of degree at most n. It is well known (see for example 
[10]) that any trigonometric polynomial / <E r„ can be uniquely determined by its values at 2n + 1 arbitrary 
points — 7r < x\ < x 2 < ■ ■ ■ < x 2n +i < 71", (we change our notation for simplicity) through the formula 

2n+l , . 
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where the polynomials t k € r„ are given by 
(17) 



2n+l , _ . 



To reach our goal we need to calculate at the nodes. The differentiation of (17) and algebraic manipu- 
lation of the result yields 

2n+ ) ( - ) 
t'(xk) Y cot h 2 ' ' 3 ' = k ' 
i=i 



t'(Xj) CSC 



3 + k, 



where 



l \ 2n+l / \ 

, , , x . {x-Xk) yj . {x-xi) 
t{x) =t k (x)sm = || sin . 



i=i 



Therefore, f'(xj) takes the form 
(18) 



' (xj - x k ) 



/'fe) - y' — — + lf(xj) y' ^t v " ' 



fc=1 sin — g — " fe=i 
Note that f'(x) is again an element of r„ and that this equation has the structure of (4) with the definitions 

D = D V = TDT-\ 



( N 



(19) 



Da = { 



EH ,{xi- xi) 
- cot . 
2 2 



1 \%i $j ') 
- CSC — — . 

2 2 



i ^ 3, 



T i3 = t'(xi)Sij, i,j = 1, AT. 



Let choose the nodes such that 
(20) 



-tt+^P J -1,2, •••,7V. 



Then, it is not difficult to see that regardless the parity of N, the product (17) and the sum given by Djj, 
evaluated at (20), satisfy the equations 



t k (x k ) t'(x k ) 



and 



N 



D jj = Y^ cot ^4^=0, 3 = 1,2,-. -,N. 



fe=i 



A' 



Therefore, substituting these expressions in (18) and taking into account(16), we obtain for this choice of 
nodes and N ~ 2n+ 1, 



(21) 



A 2lT T 
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where we have defined the TV-dimensional matrix L z as 



(22) L z = -iD v 

(we take K = 1). Therefore, the equation (14) for the discrete rotation operator A becomes 

2tt 

(23) A = e~N r v =e- zeL *, 



where e = 2tt/N. Like it comes out in the continuum case, the argument of the exponential is proportional 
to the angle of rotation (in this case e) times the derivative D v . Thus, the interpretation of L z given by (22) 
as generator of discrete rotations follows immediately. On account of (21), (13) and (15), L z has the 2n + l 
eigenvalues 

{—n, —n + 1, • • • , —1,0, 1, • • • ,n — l,n}, 
and the normalized eigenvectors |m), m = — n, ■ ■ ■ , n whose components, in the </?-basis are 



(ifk\m) = 




-im<p k 



where ip^ = fee, k = 1, 2, • • • , N. 

The fact that the quantum rule for the ^-component of the angular momentum can be projected in 
a finite-dimensional space maintaining the form it has in the continuum case (a consequence of the use of 
D v as a projection of ■£-), reinforces the possibility of the construction of a finite-dimensional algebra for 
Quantum Mechanics as it appears in other problems [11]— [12]. 

We consider now the case TV = 2n, but before we have an important remark. According to the 
uniqueness of the representation of trigonometric polynomials in terms of the Dirichlet kernel evaluated at 
differences of (20) (see for example [10]) we have that D v equals to (2n + l)/2 times the derivative of the 
Dirichlet kernel evaluated at Xj — x k . On the other hand, it should be pointed out the dependence of the 
explicit form of the matrix representation of on the type of projection functions; this means that it is 

possible to construct other matrices representing the derivative if we restrict the points to be in (0,7r] 
and project on the cosine polynomials. 



3.2. Even case 

Essentially, the case of N = 2n differs from the odd case only by the space of functions where the 
projection takes place. Let us begin by considering a function of the form 

(24) g(x) = sm(x/2)f(x), 

where / e r n _i. Let — ir < x\ < x* 2 < ■ ■ ■ < x* 2n _ x < ir be N — 1 arbitrary points, different from zero. We 
can interpolate / at these nodes to yield 

to) - E/(4) Sin( !/^ W - E [sin(4/2)/(4)]-f (X/2) ^ 



t%{xl) " K ' ^ J sin(4/2)t*(x*)- 



fc=l KV K/ k=l 



that is, 



X ' sm(x/2)t* k (x) 



( 5) 9{x) ~ ^/ Xk KHxi/2)tUxiy 

The functions t%(x) are Gaussian polynomials, given by a product like (17) at the N— 1 nodes x*. Now, let us 
consider the set of points formed by zero and Xj, j = 1, 2, • • • ,N—1, and denote them by Xj, j = 1, 2, • • • , 2n. 



Let tk(x) be the corresponding basic interpolatory polynomials. Then, taking into account that g(0) = 0, 

(25) becomes 

(26) 9(x) = | 9(xt) ^, 

and we have an interpolation formula for functions of the form (24) at the points Xj. Therefore, most of the 
reasonings made in the odd case can be applied to (26) to yield a formula similar to (18) for g(x), so that, 
for 2n points given by any 2n — 1 different points of (— tt, tt] and zero, the matrix 

(27) D = D V = TDT-\ 

with the definitions (19) is an exact representation of the derivative of functions given by (24). Moreover, 
taking into account that cosa;/2 = sin[(x + 7r)/2] and that tt can be taken as one of our nodes and that this 
function vanishes just at x — tt, we can manipulate the function h(x) = cos(x/2)/(x) in the same form as we 
did with (24) to conclude that the fc-th power of (27), constructed with a set of 2n distinct points of (—tt, tt] 
that includes zero and tt, is a 2n-dimensional matrix projection of the fc-th derivative on the subspace of 
functions of the type 

e ix/2 f(x), 

where / £ T n -i- Now, if we restrict these nodes to be equidistant as in (20), the same formulas (21)-(23) for 
the relations between A, D v , L z , and A are yielded, but now for N = 2n. This case makes up a two-valued 
or spin representation of the finite group of rotations. 

Using the same definitions of the preceding case, we have that L z has the 2n eigenvalues 

{-(2n - l)/2, -(2n - 3)/2, • • • , -1/2, 1/2, ■ • • , (2n - 3)/2, (2n - l)/2}, 

and the normalized eigenvectors \m), m — — (2n — l)/2, • • • , (2n — l)/2 are given again by 

{<p k \m) = ^=e- im *>>, 



where iph = ke, k = 1, 2, • • • , TV. 

To end this section let us remark that the validity of formulas (21)-(23) for both N odd and even shows 
one discrete formulation of rotations in terms of an finite- dimensional matrix representation of the derivative 
of certain periodic functions that can be associated, according to quantum postulates, to one component of 
the angular momentum operator projected in C N , giving thus, a finite subgroup of the rotation group that 
shares some of the properties of the full group. 

4. The numerical eigenproblem of L 2 . 

Our purpose in this section is to obtain a finite-dimensional matrix representation of the angular 
momentum eigenproblem of a system described by three classical degrees of freedom by applying the results 
of sections 2 and 3. 

As usual, we choose the spherical variables 9 and ip to describe the problem. In order to apply the 
results of section 2, we need two matrices: Dq of dimension N, and D v of dimension M, to represent d/d6 
and d/dip respectively. Since —tt < ip < tt and the functions to represent (the spherical harmonics Y™(9, ip)) 
are trigonometric polynomials, we will use as matrix D v , the one given by (19) with the change of notation 
Xj — > ipj, where the M points <pj are given by formula (20): 

Concerning the variable 9, we have several alternatives to choose a N x N matrix for Dg. Among these, we 
present only two of them yielding exact results at the nodes. The first one follows the ideas given here and 
the second is taken from [13]. 
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4.1. A matrix for L 2 

Since < 9 < ir, we can take any set of N distinct points Oj of (0, ir) to construct the matrix Dg 
according to (19). Besides, the fact that Y™(6, ip) are defined for integral indexes excludes the use of the 
spin representations, i.e., the matrices D v and Dg should be constructed with an odd number of points. 

Now, according to the results of section 2, the differential eigenvalue problem for f s (0, p) = Y™(0, ip), 

d 2 f s ^ ..df, 1 d 2 f s 

W + cot + ^Te W = sfs = ~ n[n + )/s ' 

takes the matrix form 

(28) L 2 f:-A:f:, a = i,2--.,™, 

where L 2 is the NM x NM matrix given by 

(29) L 2 = [D e 2 + cot(0)D e + sin- 2 (#)D/] , 

j* e qNM an( j ^* j g j n g enera i a com plex number. The matrices of (29) are given by 

9 = 1 M O ©, De = 1m ® ^e, D v = D v ® ljv, 

where is a diagonal matrix with entries Oj^ = OjSjk, and ljf is the identity matrix of dimension K. Thus, 

(29) can be rewritten as 

(30) -L 2 = 1 M ® [^e + cot(0)D o ] + Z? 2 ® sin- 2 (6). 

Since <p) is the (tensor) product of the associated Legendre functions P™{0) and e lm6 and, 

on the other hand, we know (from the results of section 3) that the degree of the polynomial that can be 
differentiated through (19) yielding exact values at N nodes (N odd) is at most (N — l)/2, formula (28) 
reproduces exactly the first 

eigenvalues and the corresponding unnormalized functions Y™(9, ip) evaluated at the nodes (9j, (pk), provided 
M > N. In this case L 2 will have necessarily (N + l) 2 /4 real eigenvalues given by 

\*=n{n+l), 

ordered according to 

s = n 2 +n + m+l, n = 0, 1, • • • , (N — l)/2, m = — n, — n + 1, • • • , n — 1, n. 
The eigenvector corresponding to A* , has the components f* s given by 

f* _ pm/f) \ imke 

■Irs — ^nm 1 n \v 'j jc , 

where e = 2n/M, c nm is a normalization constant, 9j G (0, it), and the relation between r, j and k, is given 
by (6). 

It is possible to choose the #-nodes in such a way that L 2 becomes a positive semidefinite matrix 
(save to a similarity transformation). To see this, note that — D 2 is positive semidefinite whereas sin _2 (0) 
is positive definite since 9j € (0, ir), j — 1, 2, • • • , N. Therefore, according to (30) we only have to find the 
conditions which make the matrix 

(31) T^LjTe = -Dg - cot(9)D e 
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positive semidefinitc. To this end, let us separate the main diagonal of Dg by writing 

Dg = TgD*T B x + d, 



where 



0, i=j, 



D* j = \ 1 {e & ) d«=^£ 2 COt 

-CSC i^j, 1=1 

[cf. (19)], and Tg is given by (19). Thus, 

(32) L 2 e = —T S [D* 2 + D*d + dD* + d 2 + cot(Q)D* + cot(6)rf]T - 1 . 

If we can find points 9j such that d~— cot(9)/2, i.e., 



N 



(33) ^' cot (^) = _ cot( ^.) 

i=i z 

[cf. (3)], equation (32) becomes 

L 2 g = Tg(—D* 2 + dD* - D*d + d 2 )^ 1 = T e (-D* + d)(D* + QTf 1 = TgDlDgTg 1 , 

where Dg is the transpose of Dg. It remains to show that a solution of (33) always exists. Since cot a; is the 
logarithmic derivative of sin a;, Eq. (33) is the condition for a critical point of the function of TV variables 

N N (zi-z) 

U(z) = U(z 1 ,z 2 ,- ■ ■ ,z N ) = Y[ sin(z fc ) J|sin - 3 , 

fe=l i>j 

where < Zj < ir, j = 1,2, ■■■,N. The existence (and uniqueness) of the solution can be proved along 
the same lines given in [14]. It is worth to be noticed that the similarity transformation is not essential 
for an eigenproblem like (28) since the similarity matrices corresponding to 9 and <p can be collected into a 
NM x NM diagonal matrix S to write (28) in the form 

Lpg s = A*g s , 

where g s = S _1 f^ and Lp is positive semidefinite. Thus, if we construct l? e with the set of nodes satisfying 
(33), L 2 is a positive semidefinite matrix (save a similarity transformation), a necessary property from the 
numerical and physical point of view in any projection scheme. 



4.2. Other matrix for L 2 

In [13], a matrix representation of the derivative for trigonometric polynomials of definite parity is 
given and, therefore, it can also be used to construct a representative of (31), yielding a matrix for L 2 with 
a higher degree of approximation than (29). The cost we have to afford for this is the lacking of simplicity: 
we can not use a single matrix for to be substituted directly in (31) unless we accept to loose precision 
in the results (see [13]). According to this scheme, the representation of a differential operator is formed 
following a given rule where certain matrix D is involved. Such a matrix is constructed with N arbitrary 
distinct points 6j € (0, ir) through a formula similar to (19) 

D = D = St>S-\ 



11 



( N 



?, = 



N 




I cot(e t -9j), 



The form that iA 



adopts in this case is 



(34) 



Lg = -D 2 , - cot(e)Do + NSOS 



-l 



where 9 is again a diagonal matrix with entries Qjk = OjSjk and O is a projection matrix with ones 
everywhere, i.e., Ojk = 1. According to [13], the degree of approximation of (34) is higher than that of (31). 
While the former yields exact results at N nodes (N odd or even) for trigonometric polynomials (of definite 
parity) of degree N, the latter produces exact results for polynomials of degree (N — l)/2 (N odd). Thus, 
the use of (34) in (30) gives the matrix 



with a higher degree of approximation: if N is an odd integer and M — 2N + 1, then (35) produces the first 



exact eigenvalues (and eigenvectors) while (30) yields only the first (N + l) 2 /4 for the same values of N and 
M. 



5. Final remarks 

Summarizing, we have found a finite-dimensional representation of the square of quantum angular 
momentum with the following properties: 

1. The coordinate representation of the ^-component of the angular momentum operator is maintained 
in this discrete scheme and the projection of L z is the generator of a finite subgroup of the group of 
rotations. 

2. The projection of L 2 is given in terms of finite-dimensional representations of the partial derivatives 
according to the well-known formula for L 2 . 

3. The spectrum of this projection contains the first eigenvalues of L 2 (counting the degeneracy) in such 
a form that the whole spectrum of L 2 can be reobtained when the number of nodes tends to infinity. 

4. The eigenvectors of this projection, corresponding to the exact eigenvalues of L 2 , can be converted 
into the exact eigenf unctions through an interpolation at the nodes. Again, this process yields the 
complete set of eigenfunctions of L 2 when the dimension of L 2 tends to infinity. 

Finally, we note that these properties makes L 2 suitable for numerical applications to quantum prob- 
lems as it will be shown in a subsequent work. 
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(35) 



-L 2 = 1 M <8> [D 2 . + cot(9)D e - NSOS- 1 } +D 2 V ® sin- 2 (9) 



N 



]T(2Z + 1) = (A+1) 2 



12 



REFERENCES 



[1] F. Calogero, Lett. Nuovo Cimento 35, 273 (1982). 

[2] F. Calogero and E. Franco, Nuovo Cimento B 89, 161 (1985). 

[3] M. Bruschi, R.G. Campos and E. Pace, Nuovo Cimento B 105, 131 (1990). 

[4] R.G. Campos, Bol. Soc. Mat. Mexicana, 3, 279 (1997) 

[5] R.G. Campos y R. Mufioz B., Rev. Mex. Ffs. 36, 1 (1990). 

[6] F. Calogero, J. Math. Phys. 34, 4704 (1993) 

[7] E.W. Cheney, Multivariate approximation theory: selected topics (SIAM Publications, Philadelphia, 

Pennsylvania, 1986) 
[8] F. Calogero, Lett. Nuovo Cimento 38, 453 (1983). 

[9] I.S. Gradshteyn and I.M. Ryzhik, Table of integrals, series, and products (Academic Press, London, 
1994), 5th. Ed. 

[10] A.F. Timan, Theory of Approximation of functions of a real variable (Dover Publications, Inc., New 
York, 1994) 

[11] R.G. Campos, Rev. Mex. Fis, 29, 217 (1983). 

[12] R.G. Campos, Rev. Mex. Fis., 32, 379 (1986). 

[13] F. Calogero, Lett. Nuovo Cimento 39, 305 (1984). 

[14] T. Popoviciu, Bull. Math. Soc. Roumainc Sci., 38, 73 (1936). 



13 



